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Abstract 

Exact solvability of the discretized N— point version of the VT— symmetric square- 
well model at all N is pointed out. Its wave functions are found proportional to the 
classical Tshebyshev polynomials Uk cos#) of a complex argument. A compact secu- 
lar equation is derived giving the real spectrum of energies at all the non-Hermiticity 
strengths Z E (—Z crit (N), Z crit (N)). In the limit Z — * the model degenerates to a 
Hermitian Hiickel Hamiltonian. 
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1 Introduction 



The essence of the concept of the so called PT — symmetric Quantum Mechanics [1] 
lies in a half-forgotten fact that although the operators of observables (say, Hamil- 
tonians H) must be Hermitian (with respect to a metric 6 in Hilbert space), they 
are also allowed to be Hermitian with respect to a nontrivial metric, 

H* = eHQ-\ / ^ = 0t > o. (1) 

In sufficient detail, the formal aspects of this idea were already well described in 
the review paper [2] showing that in nuclear physics a decisive simplification of 
certain Schrodinger equations may be achieved after a formal transition from the 
most common @i = / to another 02 7^ /. 

Beyond the area of nuclear physics the feasibility of the necessary calculations 
represents a key technical challenge for 0^1. Fortunately, for some differential 
Hamiltonians H = p 2 + V(x), an unexpectedly satisfactory answer has been found 
in a factorization of 0, typically, into a product of parity P and the so called quasi- 
parity Q [3] or charge C [4]. In the other words, one requires that the observables 
are both VT— symmetric and CP— pseudo-Hermitian. This means that all our non- 
Hermitian observables must be compatible with eq. (1) and that they must also com- 
mute with certain operator PT. In the terminology advocated by A. Mostafazadeh 
[5] the latter requirement should be generalized and re-interpreted as the so called 
P— pseudo-Hermit icity relation 

H j = PHp-\ I^p = p^ (2) 

In such a widely accepted scenario [6] our observables may remain manifestly non- 
Hermitian in the current sense, H 7^ if'. Still, their spectra must remain real and the 
work with the underlying indefinite 'pseudo-metric' operators of generalized parity 
P should remain sufficiently easy. 

In the light of these two requirements we feel particularly inspired here 

• [a] by the rigorous Krein-space analysis [7] of the reality of spectra of quantum 
particles confined inside a one-dimensional "PT— symmetric box V(x), 
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[b] in context [a] , by our older explicit construction [8] of bound states in one 
of the most elementary square- well forms of V(x), 

[c] in context [b], by the very recent Weigert's [9] three-by-three matrix dis- 
cretization 



2 + |iZ -1 







-1 2-\iZ 



7 

V ^ J 



4 



7 

V Po J 



of the oscillator of ref. 



(3) 



In this general framework we intend to start from the VT— symmetric ordinary 
differential Schrodinger equation 

d 2 



dx 2 



+ V(x) 



i})(x) = Eip(x), 



^(±1) = 0, 



V(x) = [V(-x)\* (4) 



and from the proof given in ref. [7] that the resulting spectrum of energies E = E n , 
n — 0, 1, ... is real and discrete for all the complex VT— symmetric potentials V(x) 
which are not too strong, 



< -7i 2 w 3.701. 



(5) 



In the spirit of item [b] we shall only pay attention to one of the simplest piecewise 
constant and purely imaginary antisymmetric potentials 



V(x) 



(6) 



+iZ x G (-1,0), 

-iZ x e (0,1). 

Along the lines of item [c] we shall introduce the Runge-Kutta discrete lattice of 
coordinates, 
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x = -1, 



Xk = Xk-\ + h = — 1 + kh, 



h 



TV' 



k = l,2, 



N 



and pay attention to the general discrete analogue 
ip(x k+1 ) - 2ip(x k ) + ip(x k -i) 



h 2 



-isign(x k )Zip(x k ) = E^(x k ) 



(7) 
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of our differential Schrodinger equation (4) + (6). In combination with the boundary 
conditions 

i/)(x ) = ip(x N ) = 0, 

this in fact represents a generalization of the Weigert's N = 4 eq. (3) to all the 
integers N. 



2 Models with the even N = 2n + 4 



Let us recollect that the continuous N = oo model (6) is exactly solvable [8]. The 
solvability in closed form also characterizes its modifications with periodic boundary 
conditions and/or more discontinuities [10]. In this context, the exact solvability of 
the following (N — 1)— dimensional matrix generalization 



i£-F -1 
-1 ■■■ 



-1 

-1 



V 



a r 



7 



Pn 



V ^ ) 



(8) 



-1 -if - F ) 

of the three-dimensional Weigert's model (3) with the re-scaled energy eigenvalues 
F = E h 2 — 2 and with the re-scaled strength f = Z h 2 of the non-Hermiticity would 
not be too surprising. 

In the first step of our analysis, due to the VT— symmetry of our problem we 

may set 



a k = a k + ib k , (3 k = a k - ib k = a* k , 



k — 0, 1, . . . , n 



with some real elements 7 = i/)(0), a k = Reip (x k +i) and b k = Imip (x k+ i). Once we 
recollect the definition of the classical Tshebyshev polynomials of the second kind 
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[11], 



Uk(cos6) = 



sin(£; + 1)0 
sin 9 



k — 0,1, ... , 



we immediately obtain the wave functions in closed form, 



oik = U k 



(a + ib), 



k — 0,1, . . . ,n . 



(9) 



This reduces the full tridiagonal (N — 1) x (N — 1)— dimensional matrix eq. (8) to 
the mere three matching conditions, 

' C/„_i (a + ib) 

Un (=?P) (a + ib) 
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' -1 if - F -1 ' 



V 



-1 




-F -1 
-1 -^~F -1 



(a-ib) 



= 0. 



. U n -i (a - ib) _ 

The first and the third lines may be simplified to give 



7 = £4 + i (a + ib) = U n+1 ( -^—^) (a - ib) . 



The middle line defines the product 



(a + ib) - U n ( F n ^ ) (a - ib) . 



2 J v ' V 2 
This forces us to separate the F = case as exceptional. 



(10) 



(11) 



(12) 



2.1 The existence of a nontrivial solution at F = 

Once we set, tentatively, F = 0, it is easy to deduce from eq. (11) that the parameter 
a must vanish for the even n — 0, 2, 4, . . . (and we may normalize 6 = 1) while 6 = 

and a = 1 for the odd n = 1,3,5, Thus, eq. (11) degenerates to the mere 

definition of the last element 7 of the eigenvector and we are left with the single 
secular eq. (12) which acquires the following two alternative forms, 



U n (±i£)-U n (±iZ)=0, n = 2m, 
U n (±iS) + U n (-±iS)=0, ra = 2m + l. 



(13) 



At any m = 0,1,... these conditions are satisfied identically. We may conclude 
that our tentative "guess of the energy" was correct and that F = is always the 
eigenvalue. It is remarkable that in spite of the non- Hermit icity of the Hamiltonian, 
this "robust" eigenvalue remains real at all the real couplings Z e (—00, 00). 



3 Closed secular equations for N = 2n + 4 

Whenever F ^ we may treat eq. (11) not only as the condition of vanishing of the 
imaginary part of 7, 

U n+ 1 (-^) (« + i& ) = Ur+l (-^) (« - i& ) ( 14 ) 

but also as an explicit definition of the non-vanishing left-hand-side quantity F7 in 
eq. (12). Its insertion simplifies the latter relation, 

Tn+i (-^) (« + i& ) = - T n + i (-^) (° - i6 ) ( 15 ) 

where T fc (z) denotes the k— th Tshebyshev polynomial of the first kind. 

One of the latter two relations defines the normalization vector (a, b) = (a , 60) 
while their ratio gives 

t„ (=Z±*) u M + r„ +1 (z£^) ^ (z£±ii) . (16 , 

This is our final secular equation which defines, in an implicit manner, the energies 
F as functions of the couplings £. 

An efficient numerical treatment of the latter eigenvalue problem may be based 
on the re-parametrization 
—F + i£ 

= cos (p, Rep = a, Imp = (3 (17) 

i.e., 

1 1 

- F = — cos a cosh f3, - £ = — sin a sinh /3 . (18) 
In opposite direction, the inversion of this change of variables 
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cosa = F, sinh/? = ^= \j F* + ? - 4 + V^ 2 + f 2 - 4) 2 + 16 £ 2 - 



transforms eq. (16) into the compact and transparent trigonometric secular equation 

Re sin[(w + l)y>] cos[(n+ l)y?*] = Q 
simp 

Its roots may be determined, numerically, as lying in the domain with negative (3 < 
and with a G (0, 7r/2) for the negative F < and with a G (tt/2, tc) for the positive 
F > 0. In this picture, the constant value of the coupling £ > is mapped upon 
a downwards-oriented half-oval in the a — (5 plane with a top at a — ir/4. Its two 
asymptotes a = and a = n/2 are reached in the limit /? — > — oo. 

In the new graphical representation the robust, £— independent energy level F = 
lies on the top of the half-oval while its decreasing and increasing neighbors are found 
displaced to the left and right, respectively, along the half-oval downwards. At the 
first few lowest N = 2n+4 the coordinates of these eigenvalues remain non- numerical, 

F = 0, F ± = ±v/2^e, n = 0, 
^0 = 0, F ± ,± = ±v/2 - e ± V^W, n=l 

etc. The closed form of these definitions enables us to determine the closed form of 
the respective critical values at which the spectrum ceases to be real, 

Z crit {4) =4^^5.66, n = 0, 
Z mt (6) = 9/2 = 4.50, n = 1, 

followed by the numerically calculated Z crit (8) ps 4.463 (at n = 2), Z crit (10) fx 4.461 
(at n = 3), Z crit (12) fx 4.463 (at n — 4) etc. These results do not contradict the 
expected n — ■> oo limit Z crit (oo) fx 4.475 derived in ref. [12]. 

3.1 A real-matrix re-arrangement of eq. (8) 

We may split eq. (8) in its real and imaginary parts and 'glue' them together in 
the following very natural pentadiagonal or, if you wish, block-tridiagonal eigenvalue 



6 



problem 



V 



i -F 


— 1 U 

-1 








— 1 
















-1 

-1 








-1 

-1 




-1 









-2 


—F 



a 



ai 
6i 



= 0. 



(20) 



/ v 7 ; 



This equation may be re-written in the partitioned-matrix notation, 
/ 



X 
-1 



-1 
X 



Co 



Cl 



V 



= 0. 



(21) 



V 7 / 



"-. '•• -1 

-1 X d 

2d T —F 

The obvious boldface two- by- two matrix elements degenerate, in an 'odd' anomalous 
last row and column, to an auxiliary vector d T = (1,0). 

A few comments are due. Firstly, all our wave-function components are now 
re-interpreted as proportional to the classical Tshebyshev polynomials Uk with a 
two-by-two matrix argument X. Equation (21) leads to an alternative formula for 
the eigenvectors, 



X 



l-F 



0,1,. ..,71+ 1. 



(22) 



V i -f j 

Secondly, once we introduce a complex angle a we may parametrize F = — g cos a 
and £ = g sin a using an optional, redundant parameter g. A peculiar feature of our 
matrices X is that their powers remain elementary in this representation, 



X m = g m 



cos ma -sinma 
sin ma cos ma 



This means that all the formulae containing polynomials (22) remain amazingly 
transparent. 

The existence of the explicit solutions (22) reduces eq. (21) to the two secular- 
equation constraints imposed upon the vector c n+1 . Of course, they are equivalent 
to our complex matching conditions as mentioned above. 



4 Solutions at the odd N = 2n + 3 



Weigert [9] did not notice that a "one-step easier" discretization of eq. (4) emerges 
at the odd N = 2n + 3, with 2N + 2 energy roots at n > 0. An alternative to eq. (8) 
then reads, in the same notation, 



i£ - F -1 




-1 ••• 




-1 




-1 


-1 


-1 






-1 




-1 -i£-F 



\ 
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ol. 



oil 



0. (23) 



Definition (9) of the eigenvectors remains unchanged but the matching condition is 
just one, 

/_F_i_iA / — F — iA 

(24) 



" = U n+1 \ ) (a + ib) = U n ( -^-^ ] (a - ib) 



2 J V 2 

The ratio between this equation and its Hermitian conjugate eliminates all the nor- 
malization ambiguities and leads to the odd— N counterpart of eq. (16), 



2 I 11 \ 2 / ~ \ 2 / ~ n ^ \ 2 / 
This secular equation is our final result. As an implicit definition of the iV = 2n + 3 
energy levels F = F(£) it possesses the compact non-numerical solutions at the first 
two values of n again, 

F ± = ±VT=e, n = 0, 

~ n = 1. 



(25) 



The respective elementary expressions for the critical constants 




f = 2.25, 
^ « 3.49, 



n = 0, 



n — 1 



are followed by the complex Cardano representation of the real Z crit (7) ps 3.946 at 
n = 2. At n > 2 one switches to a purely numerical algorithm giving Z crit (9) 4.148 
at n = 3 etc. In comparison with the parallel results sampled in section 3 at even N 
we notice a slowdown of the numerical convergence towards the n — > oo limit. 
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